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We discuss rotating wormholes in general relativity with a scalar field with negative kinetic energy. 
To solve the problem, we use the assumption about slow rotation. The role of a small dimensionless 
parameter plays the ratio of the linear velocity of rotation of the wormhole's throat and the velocity 
of light. We construct the rotating wormhole solution in the second order approximation with respect 
to the small parameter. The analysis shows that the asymptotical mass of the rotating wormhole is 
greater than that of the non-rotating one, and the NEC violation in the rotating wormhole spacetime 
" is weaker than that in the non-rotating one. 

^ ■ PACS numbers: 04.20.Jb, 04.25. Nx, 98.80.Cq 

(N : 

C/^ ; I. INTRODUCTION 

\ Wormholes are usually defined as topological handles in spacetime linkingwidely separated regions of a single 
■ universe, or "bridges" joining two different spacetimes 0, As is well-known they can exist only if their throats 
Q ■ contain exotic matter which possesses a negative pressure and violates the null energy condition. The search of 
O^i: realistic physical models providing the wormhole existence represents an important direction in wormhole physics. 
^ ■ Various models of this kind include scalar fields, wormhole solutions in semiclassical gravity, solutions in Brans-Dicke 
theory, wormholes on branes, wormholes supported by matter with exotic equations of state, such as phantom energy, 
the Chaplygin gas, tachyon matter, and others [1, 
\ It is worth being noticed that most of the investigations deal with static spherically symmetric wormholes because 
J> • of their simplicity and high symmetry. At the same time it would be important and interesting from the physical point 
C'l ' of view to study more wide classes of wormholes including non-static and rotating ones. Non-static wormholes whose 
I geometry is depending on time have been discussed in the literature. In 1993 Roman explored the possibility that 
inflation might provide a mechanism for the enlargement of submicroscopic, i.e., Planck scale wormholes to macroscopic 
'""J . size. He used the Hue element with the exponential scale factor. Kim Q generaHzed the Roman's consideration by 
0^ ' using the scale factor in a general form. Various aspects of non-static wormholes conformally related to static wormhole 
' geometries were investigated in [8, 9|,[lDi[lll; in particular, some issues concerning WEC violation and traversability in 
these time-dependent geometries were discussed. Kuhfittig p3| considered a spherically symmetric wormhole spacetime 
with the metric whose components are time-depending. Exact solutions describing cosmological evolution of scalar 
field wormholes were obtained in [l3. [l^. 

Rotating wormholes have also been an object for study. Some general geometrical properties of stationary rotating 
wormholes have been first analyzed by T eo fl5l| . General requirements for the stress-energy tensor necessary to generate 
' the rotating wormhole were discussed in . The WEC violation and traversability in the rotating wormhole spacetime 
were in details studied in [l3|. Kim [lB| investigated scalar perturbations in a particular model of the rotating wormhole. 
In Ref. (l9l | the authors studied a slowly rotating wormhole surrounded by a cloud of c harg ed particles. Arguments in 
favour of the possibility of existence of semiclassical rotating wormholes were given in [20|- 

In previous work [2l| we were continuing the study of rotating wormholes. Our aim was to construct an exact 
solution describing these objects in general relativity with a scalar field. As is well-known (see [H, 1131) ^ scalar ghost, 
i.e. a scalar field with negative kinetic energy can support static spherically symmetric wormholes. Moreover, such 
wormholes are stable against Hnear spherically symmetric perturbations [131 ■ In j2l| we looked for rotating wormholes 
supported by the scalar field with negative kinetic energy. To solve the problem we supposed that a wormhole is 
very slowly rotating and constructed a solution in the first order approximation with respect to a small parameter 
characterizing the velocity of rotation. In this approximation the only term ~ fldtdtp is added to the initial non-rotating 
wormhole metric, where Q is the local angular velocity of rotation. This term results in the well-know dragging effect 
in general relativity. Namely, it was shown that a test particle initially propagating along the radial direction turns 
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out to be involving into the wormhole rotation so that after passing through the throat of a wormhole it continues 
its motion along a spiral trajectory moving away from the throat. Propagation of light exhibits similar behavior. The 
ray of Hght after passing through the rotating wormhole throat is propagating along the spiral. At the same time, the 
first order approximation does not give an answer to a number of interesting and important problems concerning the 
rotating wormhole mass, the NEC violation and others. For this reason, in this paper we construct and analyze the 
second order solution describing the rotating wormhole in general relativity with the scalar field with negative kinetic 
energy. 

The paper is organized as follows. In Section|lllwe give some general formulas and write down the field equations. A 
static spherically symmetric wormhole is briefiy discussed in Section IIIII In Section IIVI we formulate the condition of 
slow rotation and introduce a small parameter characterizing the rotation velocity. Then we construct a solution 
describing a rotating wormhole in the second order approximation with respect to this small parameter. Some 
properties of the solution is analyzed in Section |Vl Namely, we discuss how a mass and a value of the NEC violation 
by a rotating wormhole differ from those of a non-rotating one. A summary of results obtained is given in Section [Vll 
In Appendix we give details of solving field equations in the second order approximation and write their solution in 
an explicit form. 



II. GENERAL FORMULAS 



Consider general relativity with a scalar field $, describing by the action 

f d^xV^ [TZ + (V$)2] , (1) 



where g^jy is a metric, g = det((;^^), TZ is the scalar curvature, and (V^)^ — g^^<^^^<^^i, is the kinetic term. Throughout 
the paper we use units G = c = 1 and the signature ( — h ++). For this signature the + sign before the kinetic term 
corresponds to negative kinetic energy, hence $ is a ghost. 

Varying the action ^ with respect to g^jy and yields Einstein equations and the equation of motion of the scalar 
field, respectively: 

■R^,, = -<^,^,<^,u, (2) 
V"V„$ = 0. (3) 

In the paper we will search for solutions of the system ([2]), ([3]) describing rotating wormholes. A spacetime with the 
stationary rotation possesses the axial symmetry. As is known (see, e.g. |2g|) a general axially symmetric metric can 
be given in the following form: 

ds^ = -Adt^ + Bdr^ + [de^ + sin^ 6{dLp - Qdtf] , (4) 

where A, B, R, il are functions of r, 9. The function Q has an expHcit physical sense; it represents an angular velocity 
of rotation in a point (r, 9) . The requirement of finiteness of the angular momentum J measured by a distant observer 
yields the following asymptotical condition for ft (2^ : 

2 / 

n^-j + 0{r-^) as r ^ C50. (5) 

Also, requiring that a spacetime should be asymptotically fiat we have A ^ 1, B ^ 1, and R^ ^ r"^ as r — > oo. 

Eqs. {SI), Q, written for the metric are second-order partial differential equations for five functions A, B, R, 
fl, and $. Solving these equations in a general form is a rather complicated mathematical problem. To simplify the 
problem, in what follows we restrict ourselves by the case of slow rotation. 



III. STATIC SPHERICALLY SYMMETRIC WORMHOLE 



To formulate the condition of slow rotation, first of all we will discuss the static spherically symmetric case. A 
static spherically symmetric solution in general relativity with a ghost scalar field was first found by Ellis |2^ and 
independently by Bronnikov jl^l. This solution can be presented as follows (see [l^): 

ds^ = -e2^('')di2 + g-2«(r)[^^2 _^ (^2 _^ a^)(^d9^ + sin^ 9d(p^)], (6) 

I a2^1/2 

^(0 - , t/2 ^ ^(^)^ (7) 
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where the radial coordinate r varies from —oo to oo, m and a are free parameters, and 

. , m / r TT\ 

u{r) — — arctan — . (8) 

a \ a 2/ 

Taking into account the following asymptotical behavior: 

r 

we can see that the spacetime with the metric ^ possesses two asymptotically flat regions. The parameter m plays a 
role of the asymptotical mass for a distant observer located at r = cxd. We will assume that to > 0. The asymptotically 
flat regions are connected by a throat whose radius corresponds to a minimum of the radius of two-dimensional sphere, 
i?^(r) — e~^"^''^(r^ + a^). The minimum of R{r) is achieved at rth — m. The value Rth — Rinh) is called the radius 
of wormhole throat. It is worth noting that there exist massless wormholes with m = 0. In this case the metric ([6]) 
takes the especially simple form: 

ds^ ^~de + dr^ + (r^ + a^){de^ + sin^ Odif^). (9) 

It is interesting that the metric ((9| was proposed a priori by Morris and Thorne in [l| as a simple example of the 
wormhole spacetime metric. 



IV. ROTATING WORMHOLE 



Now let us consider rotating wormholes. For this aim we take the wormhole metric in the form (|4]) with r £ 
(— cx),+cx)). Assume that the throat of a wormhole corresponds to the value r = rth- Deflne the throat's radius as 
-Rth = -R|r=rth.e=7r/2- Also we introduce the value Slth = ^^|r=rth,e=7r/2j being the equatorial angular velocity of rotation 
of the wormhole throat. Without loss of generality we may suppose that Jlth > 0, i.e. the throat is rotating in the 
positive direction. Assume now that the following condition is fulfllled: 

i?th^^th < c, (10) 

where c is the velocity of Hght. The condition (fTO|) means that the linear velocity of rotation of the throat is much 
less than c. Further, we will consider an approximation of slow rotation with the small dimensionless parameter 
A — i?thf^tii/c. In this approximation components of the metric |[4l), describing the rotating wormhole, should just 
slightly differ from respective components of the static metric Following the procedure given in [IBl, we represent 
the metric functions A, B, R, Q and the fleld $ as an expansion in terms of powers of A: 

n ^ Xuj + 0{X^), (11) 

A = Aoil + X^a) + OiX*), (12) 

B = Bo(1 + A2/3) + 0(A^), (13) 

R = i?o(l + AV) + 0(A4), (14) 

$ = $o(1 + A20) + O(A4), (15) 

where Aq, Bq, Rq, and $0 are zero order solutions corresponding to the unperturbed static spherically symmetric 
conflguration (tSJ, Jll): 

Ao^e'-, B, = e-^\ R^ ^ e'^^ir^ + a'), $„ ^ l!^^!±^ ^(r). (16) 

ztt'-'^ to 

It is necessary to emphasize that the substitution A —A (or, equivalently, ilth —^th) merely corresponds to 
the rotation in the opposite direction. It is obvious that in this case the angular velocity fl is also changing its sign, 
fl while the functions A, B, R, and $ do not depend on the direction of rotation. Mathematically, this means 

that f2 is an odd function of A, i.e. il(— A) = — f2(A), and the others are even, i.e. ^(— A) = ^(A), etc. Therefore, the 
expansion ifTTj) for Q contains only odd powers of A, while the other expansions (|12m5p contain only even powers of 
A. 

Substituting the expansions I|11H15P into the fleld equations (121); ([SI and collecting terms with similar powers of 
A yields a sequence of n-th order equations with n corresponding to powers of A. The zeroth order fleld equations 
describe the static spherically symmetric conflguration Cutting the sequence of equations on a deflnite n 

corresponds to the n-th order approximation of the theory. 
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A. The first order approximation 

Rotating wormholes in the first order approximation have been studied in [2l|. As is seen from Eqs. (fTT1) - lfT5|) . in 
this approximation the functions A, B, i?, $ remain to be unperturbed, while the angular velocity fi, being initially 
equal to zero, takes the form Q ~ Xoj. The only nontrivial equation for uj reads 

^ de[sm^ e dgoj] = (r^ + a^)d^^uj + 4(r - m)drUj. (17) 



sin^ u 

Its solution, having a natural physical meaning, is 

uj{r) = 

a 

with 



(18) 



LOo{^i) = [l-e-'^^^{l + 8^i')]-\ (19) 

where = m/a is a dimensionless mass parameter. This solution describes a slowly rotating wormhole with the 
angular velocity fl = Xuj and the angular momenta J± : 

J± = l^iu;oq±{a^ +4m^), g+ = 1, g_ = g-^-f , (20) 

which are defined from asymptotics uj — 2J±\r\^^ + 0(|r|^^) at r ^ zLoo. 

Stress that the first order approximation lets only to determine the angular velocity Q of the wormhole's rotation 
and, as a consequence, reveal interesting features in a motion of test particles and a propagation of light in the rotating 
wormhole spacetime [2l||. At the same time, since the metric functions A, B and R and the scalar field $ remain to 
be unperturbed in the first order approximation one cannot answer a number of important questions. For example, 
how does the rotation changes such characteristics of the static wormhole as its mass and throat's radius? How does 
the value of violation of the null energy condition depend on the wormhole's rotation? To answer these questions, 
hereafter we will consider the second order approximation. 

B. The second order approximation 

Substituting the expansions lfTT|) - lfT5|) into the Einstein equations ^ and the scalar field equation ^ and collecting 
A^-terms yields the following system of equations for a, /?, p, and 4>: 

[r^ + a'^)dla + {d^a + cot 9 deo) + (2r + m)dra - m{dr(3 - idrp) = e-'^''oj'^{r^ + a^f sin^ 9, (21) 

(r^ + a^){dla + Ad^p) + (9^/3 + cot 9 deP) + Sto^..^ + (2r - m)(49^p - 9^/3) 

= e-4«,^'2(^2 + a^f sin^ 9 + + g^^^^)^ ^22) 

m 

(r^ + a^){dlga + 2d^gp) - r{dea + de(i) + 2mdea = + ^g^^^ ^23) 

TO 

2{r^ + a^)dlp + 2{dlp + cot 9 dep) + d^a + d^P + (r - m){dra - drP) + 4(2r - m)drp - 2/3 + 4p = 0, (24) 

2{r^ + a^)dlp + 2{dlp + cot6'aep) + cot 61 (^ea + ^e/?) + (r - m){dra - dr(3) + 4(2r - m)drP 

-2/3 + 4p^ -e-4"tj'2(r2 + a^f sin^ 9, (25) 

2u{r^ + a^)dl(t> + 2u(a|?!) + cot 9 decf) + 4(ru + m)^^?!) + TO(ara - dr(3 + 49^^) = 0, (26) 

where drOt = da/dr, &f.a — d^a/dr^, etc. One may simplify the system l(22 |) - ((26)) noting that Eq. ([23]) can be 
integrated straightforwardly resulting in 

(r^ + a^){dra + 2drp) - r{a + /3) + 2TOa = ^ ' — '-ucj) + /i, (27) 
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where fi{r) is an arbitrary function of r. Also, it will be convenient to consider their combinations [ l(24|) + l(25|) ] and 
[{SI) -{Ml)] instead of equations ^ and jM]): 

4(r2 + a'^)d^p + A{d^p + cote dep) + (d^a + cot Odea) + {d^(3 + cotOdeP) + 2(r - m){dra - dr(3) (28) 
+8(2r - m)drp - 4/3 + 8p -e-^^'iv'^ir^ + a^f &\t? 6, 

dla + dlP - 0016* (aea + dep) = e^'^"w'2(r2 + a^f sin^ 9. (29) 
Then, integrating Eq. l(29|) yields 

a + f3^ ie-4«^'2(r2+a2)2(2cos26i-l) + /3Cos6' + /2, (30) 

where /2(?') and /3(r) are arbitrary functions of r. Taking into account that the rotating wormhole configuration 
should possess the symmetry 9 ^ n - 6 one should set /3(r) = 0. Finally, the equations ([21]), ([221, |[26l, ll27|, l(28|) 
and (|30| form the system to be solved. Note that only four equations of this system are independent since the Bianchi 
identity V^G(^ = and the conservation law V^T^ — take place. 

To solve the system of field equations we will follow Ref. |25j] and expand the functions a, (3, p, and tj) in spherical 
harmonics:^ 

a(r,6l) ^ao{r)+a2{r)P2{e) + ..., (31a) 
(3{r,e) ^ fJoir) + (32(r)P2{0) + (31b) 
p(r, 9) = po(r-) + P2{r)P2{0) + (31c) 
(j){r,9) ^ Mr) + Mr)P2{0) + .... (31d) 

We would like to stress that the expansions l|3ip contain only even spherical harmonics being symmetric with respect 
to the equatorial plane 9 — ^. Another convenient simplification of the metric may be made here. Transformations 
of the type r ^ f{r) do not change the form of the metric (0]). Such a coordinate transformation may, therefore, be 
used to provide the additional condition 

Po{r) = 0. (32) 

This will be assumed in the following. Now substituting Eqs. ^ into dST]), ([221), ll27|), {SH]) and dSO]) we find 

(r^ + a2)c,^ + (2r + m)a',^ - m{(3'^ - 4p^J - n{n + l)a„ = f ^.-^"(r^ + a'fuj'^iSno - Sra), (33) 

(r^ + a'){a';, + 4^:1) + Sma^ - (2r - m){P'^ - Ap'J - n{n + l)/3„ = |e-^"(r2 + a^fu'^iS^o - Sra) 

+ , , , , , 

+8 U0n ', 34 

m 

{r^ +a^){a'^ + 2p'„)- {r-2m)an-r(3n = 4 u0„ + ./i<5„o, (35) 

m 

4(r2 + a^)p'; + 8(2r - m)p',^ + 2{r - m){a'^ - + 8p„ - 4/3„ 

-n(n + l)(4p„ + «„+/?„) = -^e-4"(r2 + a2)V2(5„o-^n2), (36) 

an+f3n - ie-4"(r2+a2)V2,5„2 + /2<5„o, (37) 
2u{r'^ + a^)(t>'^+4{ur + m)(t)'^~2n{n + l)u<j)n + m{an- I3n + 'ipny = 0, (38) 



^ Spherical harmonics or Legendre polynomials Pn{9) obey the equation (see, for example, [27| ) 



They satisfy the Rodrigues' formula 



cote — - = -n{n + l)Pn. 



P^(g) = LJ)1 ( " {[1 - cos2 0]"}. 

^ ' 2"n! \dcoseJ ' ' 



In particular, Po{d) = 1, Pi{0) = cos9, P2{6) = | cos^ ^ ~ |- Note that P„{8) = {—l)"Pn{n — 0), i.e. even spherical harmonics (with 
even numbers n) are symmetric with respect to the equatorial plane ^ = , while odd ones are antisymmetric. 
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where a prime means the derivative with respect to r, n = 0, 2, and Snii is the Kronecker delta. Note that only the 
n = and n = 2 equations involve the angular velocity ui. The coefficients in the expansion of a, (3, p, and <j) with 
n > 4 must, therefore, vanish since they vanish when the wormhole is not rotating, i.e. 

an = Pn = Pn = 0n =0, n> 4. (39) 

This reduction in the number of values of n from infinity to 2 is the central simplification of the slow rotation 
approximation. In place of a system of partial differential equations one now only has ordinary differential equations 
for the seven unknown functions ag, Po, (j^o, and a2, P2, P2, 4>2- A solution of these equations in an explicit analytic 
form is given in the appendix. Generally speaking, this solution depends on several constants of integrations. However, 
their values are fixed if one assumes that the perturbations are everywhere regular and obey the natural boundary 
conditions: 

ao|r^±oo = 0, /3o|r-*±oo = 0, (/'o|r-»±oo = COUst, 
a2|r^±oo=0, /?2|r^±oo=0, P2|r^±oo=0, 02 1 r^±oo = 0. (40) 

These conditions guarantee that there is no rotation far from the wormhole throat. Note that (t>o{r) tends to a constant 
as r ^ ±oo because the action ([1]) is invariant with respect to the shift $ ^ $ + const. In Figs. [1] and [2] we give 
graphical representation for ao, Po, 4'q, and a2, P2, P2, </'2- 




-20 -10 10 20 




Phc. 1: The graphs of ao{r), (3o{r), and (/f>o(r) for a = 1. The thick, middle and thin curves correspond to m = 0, 0.5, 1, 
respectively. 



V. ANALYSIS OF THE SOLUTION 

Properties of the non-rotating wormhole spacetime with the metric ([6]) are only determined by two parameters m 
and a, where m represents the wormhole mass measured by a distant observer located at r = 00, and a determines the 
radius of the wormhole throat. As was shown in the previous section, in addition to m and a the rotating wormhole 
solution depends on the parameter A = Rth^th/c being a dimensionless linear velocity of rotation of the wormhole 
throat. In this section we will discuss the problem: How do characteristics of the rotating wormhole differ from those 
of the non-rotating one possessing the same parameters m and a? 
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Phc. 2: The graphs of a2{r), P2(t), P2{r) and 4'2(t) for a = 1. The thick, middle and thin curves correspond to m = 0, 0.5, 1, 
respectively. 



A. The mass of rotating wormhole 



To find the rotating wormhole mass M, we should consider the limit gtt\r^oo ^ — (1 — 2M/r). In the second order 
approximation we have found gu — A{r,9) = —e'^^^^^\\ + }?[aQ{r) + a2{r)P2{9)]}. Note that a2{r)\r^oc ~ hence 
the n = 2 solution does not give any contribution in M . Taking into account the following asymptotical property: 



ao(f')|r 



2Ar 



0{r- 



(41) 



with 



2(1 + 10^i^){3^n^ + (u(/x) + ^){e-'^^''' ~ 1)) 2^le^<^'^34^I^ - A*^ + 1) 



37r 



3(1 + ^2) 



16/Lt' 



where fi = m/a is a dimensionless mass parameter, we obtain 

M = m + A^Ato(^), 



(42) 



(43) 



The value of A^Am(/Lt) characterizes the difference between the rotating and non-rotating wormhole masses M and 
TO, respectively. In Fig. [3] the graph of Am{fi) versus fi is shown. Note that Ato(^) is positive for all fi and is the 
greater the greater /i. It is worth to stress that Ato(/z) is not equal to zero in case /l( = 0: 



Am(0) 



37r 



(44) 



The case ^ — or, equivalently, m = corresponds to the massless non-rotating wormhole. The non-zero value Ato(O) 
means that there do not exist massless rotating wormholes. 



B. The NEC violation 



A violation of the null energy condition (NEC) in the vicinity of the wormhole throat is an essential feature of 
wormhole physics 0,0, Si- The NEC reads T^^fc^fc" > 0, where T^^ is the stress-energy tensor, and is a null vector. 



8 




Phc. 3: The graph of Am(/^) for a = 1. 



Using the Einstein equations, it can be represented in the geometrical form R^^k^^k^ > 0, where is the Ricci tensor. 
To analyze the NEC in the rotating wormhole spacetime with the metric Q, we choose k'^ = {A~^/^ , B~^^^ ,0, QA~^^'^) 
and introduce the value S = R^^k^k'^ ■ In the second order approximation S takes the following form 



where 



and 



eo(r) 
6M 



=2«(r) 



(r2 + a?Y 



=2n(r) 



S(r,0) = ^o{r) + X%o{r) + Ur)P2m. (45) 



SoW = -2e-W/^^, (46) 



2/3o(to^ + a^) + {r- m)(r2 + a2)(c,^ + p'^) ^ (47) 
2/32(to' + a') + (r^ + a2)[_3(a2 _ fi^) + (r - m)(a'2 + /^j - M) - 2p^'(r2 + a^)]! . (43) 



The value So characterizes the configuration without rotation. As is seen, So is everywhere negative, hence the NEC 
is violated in the whole spacetime of the non-rotating wormhole. As for rotating wormholes, it will be convenient to 
average the quantity S(r, 9) over all directions: 



S(r) 



47r7o Jo 
So(r) 



S(r, 9) sin 9d9dip 



2/3o(m2 + a^) + (r - m){r^ + a^){a'f^ + /3^) 



(49) 



In Fig. ^ the quantities S(r) and So(r) are shown together. It is seen that the value of S(r) is everywhere negative 
but greater than So(r). This means that the NEC violation in the rotating wormhole spacetime is weaker than that 
in the non-rotating one. 



VI. SUMMARY 



We have constructed a solution describing slow rotating wormholes in general relativity with the scalar field 
possessing negative kinetic energy. The role of a small dimensionless parameter A plays the ratio of the linear velocity 
of rotation of the wormhole's throat and the velocity of light, A — Rth^th/c. The field equations have been solved 
in the second order approximation with respect to A. It is worth noting that we succeeded in finding a solution in 
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Phc. 4: The graphs of H(r) (thick line) and Eo(r) (thin line) for a = 1. 



an explicit analytical form. Its analysis has shown that a mass of a rotating wormhole is greater than that of a 
non-rotating one. As a consequence, this means that rotating wormholes, in contrast to the non-rotating ones, cannot 
possess a zero mass. The respective analysis of the NEC violation in a rotating wormhole spacetime reveals the fact 
that it is slightly weaker than that in a non-rotating one. 



Acknowledgments 

The work was supported by the Russian Foundation for Basic Research grants No 08-02-00325, 08-02-91307. 



Appendix 

Here we will solve the system of equations l(33|) - (|38l) . Note that the equations for different values of n are not 
coupled together. For this reason, we will consider the equations for n = and n = 2 separately. 

The case n = 0. In this case the equations l(35l) and l(37|) turn into identities because fi{r) and f2{r) are arbitrary 
functions. As a result, we obtain 



{r^ +a^)a'' + 3ma'a-{2r-m)P' = ^e-^-^{r^ + a^fco'^ + 8 (51) 

(r-m)(ao~/3o)'~2/3o = -^e-''^{r^ + a^co'^ (52) 
2w(r2 + a2)(/)o + 4(ur + m)^;, + m(ao - /3o)' = 0. (53) 

Note that the first order equation (|52l) plays the role of a differential constraint. A general solution to the system 
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(|50l) - l(53|) can be given explicitly as follows: 



ao{x) 



(j)o{x) 



6{x — fi) 

8g4«(a;) 



1 



6/i(x — /i)-^ 



C'l + C2X + C'ixu{x) 

(1 + 10/i2)a;3 + 2^x2 + (1 + + i6/)a; - 4^3(1 + 4^2^ 
- (1 + xfi){C3fiu + + C'l) + C^ii^ip - x) 



u(x) 



3Ai(l + 2/x2)x3 - (1 + 2/^2 - 8/)x2 + 3Ai(l - 2/i2)x - - 28/ -2^-1 
1 



^lCi + m(.x)C5 + 



Ci[l + (a:; — /i) arctanx] + C2/i + C3iiu{x) 



3(2; -/(a;2 + 1) 



12(a;- ^) 

1 + 16/ - 2^(1 - 8/)a; + (1 + 4/)x2 - 2^a;^ 



(54) 



(55) 



(56) 



where x = r/a is a dimensionless radial coordinate, fj, — m/a is a dimensionless wormhole parameter, u{x) — 
^(arctana; — n/2), and Ck {k — 1, ...,5) are constants of integration. Stress that values of Ck are not free. They are 
connected by means of the constraint (|52ll . Also, the constants Ck should be chosen so that to provide an appropriate 
asymptotical behavior of the solutions: 



ao|r-»±oo — 0, (3o\r^±oo — 0, 

0o|r— ±00 = const. 



(57) 
(58) 



Moreover, a choice of Cfe should guarantee regularity of the solutions ao{x), Pq{x), (j>Q{x) at the point x — fi. Altogether, 
these conditions let us fix all values of the constants Ck as follows 



Ci = 8At(l + 10/) 



1 _ <H1 _ 1) 

TT/i ^ ' 



8Aie4"('')(34/-/ + l) 
/Tl 



C2 - -8(1 + 10/), C3 = —(1 + 10/)(e-4'^^ - 1), 

32/ -^Ci ^ C3/-Ci(l + /) 
^4 = ;r:; , O5 = 



24/i 



12/i(l + /) 



The case n = 2. The system l(33 |) - (|38| now takes the following form: 

(r^ + a^)a'i + (2r + m)a'2 - m(/3^ - 4/) - 6a2 
(7-2 + a2)(«^' + 4p^') + 3,„c.^ _ (2r - m)(/3^ - 4/) - 6/32 

(r^ + a2)(a;, + 2^^ - (r - 2m)a2 - rp2 = 4 



-2w + 8 



m 

U<p2, 



{U(f>2 



A{r^ + a^)/ + 8(2r - m)/ + 2(r - m)(a2 - /^a) - 16p2 - 602 - 10/32 = 2u, 

^2 + /32 = 

2u(r2 +0^)0^' +4(Mr + m)(?!)^ - 12M02 + m(a2 - /32 + 4p2)' 0, 



(59) 

(60) 
(61) 

(62) 

(63) 
(64) 
(65) 
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where v 



a2{x) 



i„-4ti^j,2 _|_ 0,2^2^/2 ^ general solution of the system ([BOll-lIMl) is 



+ Ufix^ + x'*(7 + 22^2) + 2Afix^{n'^ + l)+x^{5 + 36^^ + 16/) 
arctanx ^iiD2{9x^ + 1) — D4{3x^ + 1) — 



3(.x2 + 1)2 
+4:fix{3 + 4/^2 - 8fi^) ~ 64/ + 14/ + 1 
fixD2{8x'^ + 7) 



+fix^Di + 



2(^2 + 1) 



L»3(3a;2 + 1) - 3xDi 



(66) 



/32(x) = ujI 



Sx"" + 12 fix'' + a;*(7 + 22/) + 24^a;^(/ + 1) + x^5 + 36/ + 16/) 



3(a;2 + 1)2 

+4^a;(3 + 4/ - 8/) - 64/ + 14/ + 1 + 16/^(1 + 4/)^ 



arctan x 



fiD2i9x^ + 1) - Di{3x^ + 1) 



+fix^Di - 3D4X - 3D3X^ -D3 + lfixD2 + l^f!LJh. 

x^ + 1 



(67) 



P2{x) = UjI 



2g4ji(a;) 



3x'= + l^nx" + x^{7 + 46/) + 36Afa;^(2/ + 1) + x^(5 + 76/ + 80/) 



3(x2 + 1)2 

+2^a;(9 + 40/ + 16/) - 64// + 32/ + 22/i2 + 1 



■ arctan x 



-^xDiipix + I) + \D, 



D2{9fix^ + 6x + n)- 2Di{3x^ + 1) 
2x{x — fi) 



2"^ -^Vf*-^ ' -^J > 4 
4^2g4u(x) 



9fix + 4 - 



x2 + 1 



^D3{3x^ + 1) - |xD4|, 



(68) 



3(x2 + l)2u(x) 

/i arctan x 



4(/ + l)u(a::) 



4(/ + l)u(a;) 



+ Afix^ + 2x^(1 + 4/) + 8tix{2fi^ + 1) + 32/ + 16/*^ + 1 
3a;^(3/l?2 - D2 - 2fiD4) + /Da - 3132 - 2^D4 
x^{2n^D4 - Di- &iiDi) + a;(-3i:'2 + 9/L»2 - &^^D4) ~ 2fiD: 



2x132 (/ + 1) 



X2 + 1 



(69) 



where Dk^ {k = 1, ...,4) are constants of integration. Taking into account the boundary conditions 

a2|r^±oo=0, /32|r— ±00 = 0, P2|r^±oo=0, 02|r^±oo=O, (70) 

and using the two differential constraints following from the system (|60| - i[65|l we can fix the values of Dk as follows 

8n 



Di = -4/i(e 



-^''^ + 1), i32 = — (e-^"''~l), 

TT 



= ^(1 - 2/)(l + e-^-A"), D4^^{1- 3/)(l + e-*-^). 
3 A-K 



(71) 
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